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The description of interacting many-electron systems in external magnetic fields is considered in 
the framework of the optimized effective potential method extended to current-spin-density func- 
tional theory. As a case study, a two-dimensional quantum dot in external magnetic fields is inves- 
tigated. Excellent agreement with quantum Monte Carlo results is obtained when self-interaction 
corrected correlation energies from the standard local spin-density approximation are added to 
exact-exchange results. Full self-consistency within the complete current-spin-density-functional 
framework is found to be of minor importance. 
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I. INTRODUCTION 

Since its introduction in 1964, density-functional the- 
ory (DFT) 1 ' 2 has become a standard tool to calculate the 
electronic structure of atoms, molecules, and solids from 
first principles. Early on, the original DFT formulation 
has been extended to the case of spin-polarized systems 3 
which also provides a description of many-electron sys- 
tems in an external magnetic field. However, in this spin- 
DFT (SDFT) framework the magnetic field only cou- 
ples to the spin but not to the orbital degrees of free- 
dom, i.e., the coupling of the electronic momenta to the 
vector potential associated with the external magnetic 
field is not taken into account. A proper treatment of 
this coupling requires extension to current-spin-dcnsity- 
functional theory (CSDFT) 4,5 in terms of three basic 
variables: the electron density n(r), the spin magneti- 
zation density m(r), and the paramagnetic current den- 
sity j p (r) . These densities are conjugate variables to the 
electrostatic potential, the magnetic field, and the vector 
potential, respectively. 

In order to be applicable in practice, DFT of any flavor 
requires an approximation to the exchange-correlation 
(xc) energy functional. The use of the local-vorticity ap- 
proximation, 4,5 which is an extension of the local spin- 
density approximation (LSDA), is problematic in CS- 
DFT: the xc energy per particle of a uniform electron 
gas exhibits derivative discontinuities whenever a Landau 
level is depopulated in an increasing external magnetic 
field. This leads to discontinuities in the corresponding 
xc potential. 6 These discontinuities then incorrectly ap- 
pear when the local values of the inhomogeneous density 
and vorticity coincide with the corresponding values of 
the homogeneous electron gas. A popular way to circum- 
vent this problem is to use functionals which interpolate 
between the limits of weak and high magnetic fields. 7 ' 8 

Explicitly orbital-dependent functionals, which are 
successfully used in DFT and collincar SDFT, 9,10 arc 



natural candidates to approximate the xc energy in CS- 
DFT for two reasons: first, they are constructed without 
recourse to the model of the uniform electron gas and 
second, they are ideally suited to describe orbital effects 
such as the filling of Landau levels. In this way, the 
problem inherent in any uniform-gas-derived functional 
for CSDFT is avoided in a natural way. 

The use of orbital functionals requires the so-called op- 
timized effective potential (OEP) method 11 to calculate 
the effective potentials. The OEP formalism has been 
recently generalized to non-collinear SDFT 12 as well as 
to CSDFT. 13 In addition, a larger set of basic densities 
has been considered in order to include the spin-orbit 
coupling. 14,15 Recent applications of the OEP method 
for atoms 13 and periodic systems 16 have indicated that 
the difference between exact-exchange calculations car- 
ried out fully self-consistcntly within CSDFT or SDFT, 
respectively, is only minor. These works have also indi- 
cated that the inclusion of correlation energies is of par- 
ticular importance when dealing with current-carrying 
states. 

In this work we consider the OEP formalism within 
CSDFT in the presence of an external magnetic field. 
In particular, we focus our attention on two-dimensional 
semiconductor quantum dots (QDs) 17 exposed to uni- 
form and constant external magnetic fields. In addition 
to the various applications in the field of semiconductor 
nanotechnology, QDs arc also challenging test cases for 
computational many-electron methods due to the rela- 
tively large correlation effects. Moreover, the role of the 
current induced by the external magnetic field is partic- 
ularly relevant in QDs 18 making them a reference system 
in CDFT since its early devclopcmcnts. 19 Therefore, it is 
interesting to examine whether the self-consistent solu- 
tion of CSDFT differs from the result obtained by adding 
the external vector potential to the SDFT scheme, which 
amounts to neglecting the xc vector potential of CSDFT. 

As expected, we find that the bare exact-exchange 
(EXX) result is not sufficient to obtain total energies 
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in agreement with numerically accurate quantum Monte 
Carlo (QMC) results, although a considerable improve- 
ment to the Hartrce-Fock result is found. However, in- 
cluding the self-interaction corrected LSDA correlation 
energies to the EXX solution leads to total energies that 
agree very well with QMC results. In addition, within the 
given approximations, our results confirm that the role of 
self-consistent calculations in the framework of CSDFT 
is only minor. In particular, we observe that accurate to- 
tal energies and densities can also be obtained by simply 
modifying the SDFT scheme by including the coupling to 
the external vector potential. Indeed, this procedure has 
been employed in the past to partially remedy the lack of 
good approximate current-dependent functional. Here, 
a validation is provided in the more general context of 
the OEP framework. 

This paper is organized as follows. In Sec. II A we re- 
view the OEP method in CSDFT. The formalism is then 
adapted to the case of QDs in magnetic fields in Sec. II B. 
In Sec. Ill A we discuss details of the numerical proce- 
dure before presenting the results of our calculations in 
Sec. IIIB. A brief summary is given in Sec. IV. 



II. OPTIMIZED EFFECTIVE POTENTIAL 
METHOD IN CSDFT 

A. General formalism 

The Kohn-Sham (KS) equation in CSDFT reads (Har- 
tree atomic units are used throughout unless stated oth- 
erwise) 



and 



\ M V + ^ A s( r ) 



(1) 



The three KS potentials are given by 
v s (r) = i, (r) + Vh (t) + v xc (r) + ± [A^r) - A*(r)] 



and 



B,(r) =B (r)+B xc (r) 



A s (r) = A (r) + A xc (r), 



(2) 
(3) 

(4) 



where the xc potentials are functional derivatives of the 
xc energy E xc with respect to the corresponding densi- 
ties, 



v xc (?) = 



B xc (r) 



8E xc [n,n\,] % 
Sn(r) 



6E xc [n,m,j p ] 
<5m(r) 



(5) 



(6) 



-A xc (r) 



SE. 



xc [n,m,j p \ 



(7) 



respectively. The self-consistency cycle is closed by cal- 
culating the density 



n(r) = ^$t( r )$ fc (r), 

k=l 

the magnetization density 

occ 

m(r) = -MfiXi^W^W ' 



k=l 



and the paramagnetic current density 

1 occ 

jfW = ^EH' (r)v$fe(r)_(v ^ (r))<I>fc(r 



fc=l 



(8) 



(9) 



(10) 



The ground-state total energy of the interacting system 
can then be computed from 

E[n,m,j p ] = T s [n, m,jj + U[n] + E xc [n, m,j p ] 

dr n(r)i;o(r) — / rfr m(r)Bo(r) 



(11) 



dr jp(r)A (r) 
^2 / drn(r)Al(r), 



where T s and U are the kinetic energy of the KS system 
and the Hartree energy, respectively. 

Gauge invariancc of the energy functional implies that 
E xc depends on the current only through the vorticity, 



= V x (jp(r)Mr)), 



(12) 



i.e., E xc [n, j p , m] = E xc [n, u, m]. 5 This immediately 
leads to the following relation for the xc vector poten- 
tial 



V(n(r)A K (r)) = 



(13) 



If one uses an approximate E xc which is given explic- 
itly in terms of the densities, the calculation of the cor- 
responding xc potentials via Eqs. (5)-(7) is straightfor- 
ward. Here, however, we deal with approximations to the 
xc energy which are explicit functionals of the KS spinor 
orbitals These functionals are, via the Hohcnbcrg- 
Kohn theorem, implicit functionals of the densities. In 
the spirit of the original OEP formalism, the correspond- 
ing integral equations for the xc potentials can be de- 
rived 13 by requiring that the effective fields minimize the 
value of the ground-state total energy (11). Therefore, 
the functional derivatives of the total energy with respect 
to the three KS potentials are required to vanish. This 
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procedure leads to three OEP equations which are most 
conveniently written as 13 



fe=l 



-Hb ^I(r)cr* fc (r) + h.c. = , 



(14) 



(15) 



fc=i 



and 



-^[$t(r)V^(r)-(v^(r)) * fe (r 

fc=i 

where we have defined the so-called "orbital shifts" 9,20 

(17) 



- h.c. = , 
(16) 



with 



Dt.= / dr'^ xc (r')$t( r ')$ fe (r') 
1 



A ac (r) $Kr')V'^(r') - ( V'*t (r') ) $ fc (r') 



^ B B 1BC (r')$t(r') ^ Jk (r')-$](r 



SE^ 



5& k (r>) 



(18) 



The orbital shifts have the structure of a first-order 
shift from the unperturbed orbital $fc under a perturba- 
tion whose matrix elements are given by d\- . Physically, 
the OEP equations (14)-(16) then imply that the densi- 
ties do not change under this perturbation. If A xc is set 
to zero, Eqs. (14) and (15) reduce exactly to the OEP 
equations of non-collinear SDFT. 12 

Eqs. (14) - (16) form a set of coupled integral equations 
for the three unknown xc potentials, and they can be 
solved by a direct computation of the orbital shifts. 12,20 
Alternatively, one can employ the Kriegcr-Li-Iafrate 
(KLI) approach as a simplifying approximation 21,22 
which is known to yield potentials which are very close 
to the full OEP ones in SDFT. In the following we utilize 
the KLI approximation in the description of a quasi-two- 
dimensional semiconductor QD 17 in an external magnetic 
field. 



B. Application to quantum dots 

The QD is described as a many-electron system re- 
stricted to the xy plane and confined in that plane 
by an external parabolic potential vq = |m*Wgr 2 with 
r 2 = x 2 + y 2 . Following the most common experimen- 
tal setup, 17 the external magnetic field is defined to be 
uniform and perpendicular to the xy plane, i.e., Bg(r) = 



V x A (r) = B e z with the gauge A (r) = B reg/2. We 
apply the effective-mass approximation with the material 
parameters for GaAs, i.e., the effective mass m* — 0.067, 
the dielectric constant, e* = 12.4, and the effective gyro- 
magnetic ratio g* = —0.44. 

In QDs the magnetization is parallel to the exter- 
nal field, i.e., these systems show collinear magnetism. 
Therefore, the KS magnetic field B s and the magnetiza- 
tion density have only non-vanishing z-components. The 
Pauli-type KS equation becomes diagonal in spin space 
and can be decoupled into two separate equations for the 
spin-up and spin-down orbitals </?fc CT (r). We further as- 
sume that the xc potentials preserve the cylindrical sym- 
metry of the problem, i.e., 



r(r) = v xc (r) ± (J,Bg*B xcz (r), (19) 



where the upper signs are for spin-up and lower signs 
for spin-down electrons, and A IC (r) = A xc {r)eg. Due 
to the cylindrical symmetry we can separate the wave 
function into radial and angular parts as ipji a (r) = 
cyLY>(H9)Rji a (r) , where the radial wave functions Rji a (r) 
arc real- valued cigenfunctions of the Hamiltonian 



H s ] n — — 



1 (&_ 

2m* \ dr 2 



l_d__f 

r dr r* 



I 

— LO r 

2 



2 

m — r 
2 



I A xc {r) 



m*c r 



± ii B m*g*B a + v H (r) + v xc<J {r)(20) 



with the total confinement f2 = \/ lUq + w 2 /4, and the cy- 
clotron frequency tu c = B§lm*c. The radial wave func- 
tions are expanded in the basis of eigenfunctions of the 
corresponding non-interacting problem, i.e., the eigen- 
functions of the Hamiltonian (20) with the Hartree and 
all xc potentials set to zero. 

As a consequence of the cylindrical symmetry, the den- 
sities are independent of the angle and thus given solely 
in terms of r = |r|. Also, only the ^-component of the 
paramagnetic current density, as the conjugate variable 
to the vector field in this direction, plays a role, i.e., 
j p (r) = (j pT (r) +j P i(r))e e . Instead of using the density 
and the ^-component of the magnetization, one employs 
the spin-up and spin-down densities. Hence, the three 
densities to be determined are n-f(r), «x( r )> anc ^ Jp( r )- 

Consequently, the OEP-KLI equations are given as a 
3x3 matrix equation which reads 



V(r)V xc {r)=K{r), 
where the potential vector is given by 

V xc (r) = (v xc j(r),v xcL (r), -A xc {r) 



The matrix T> reads 



V 



n T (r) 

ni (r) 
jpl( r ) 3pl(r) 




(21) 



(22) 



(23) 
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where the densities and current densities are given by 

(24) 



]p<7 



{30 



The last component N{r) in Eq. (21) reads 

occ , 2 



(25) 



(26) 



The right-hand-side of Eq. (21) contains functional 
derivatives of the xc energy. They can be calculated once 
an approximation to the xc energy is specified. Here, we 
use the EXX approximation to E xc , i.e., 



E 



EXX 



-. occ 

\ E E 



o=T.l {jl},{km} 

;2 / ^ CT ( r ')^; CT (r)^ m( T(r')^fe ma (r) 

drdr — ; ; 

e* r - r' 



(27) 



The first two components of TZ on the RHS of Eq. (21) 
are then given by 

occ 

^1,2(0 = -- E R]ic{ r )Rk m <j{r) 

{jl},{km} 

x [ d 2 r > eie ' {l ~ m)R Mr')Rkm<r(r') 



\/r 2 + ?' /2 — 2rr' cos 0' 

- occ 

- -E^WW+c.c., (28) 

{3'0 

where for 7?.i cr =t and for 7?. 2 f =4- The third compo- 
nent is given by 



/ + m 



Rji a (r)R kma (r) 



<y=hl{ji},{km} 



x / d 



Vr 2 + t*' 2 — 2rr' cos ^ 

2 occ 

o E E^W%+ Ra ( 2Q ) 

°-=T4 {iO 



with 



= Jd 2 r L(r) + -^«)4(r) 



E / / d 2 rd 2 r'e ie '('-" 1 ) 

7TUT \ I ) 



\fr 2 + 7* /2 — 2rr 7 cos ^ 



(30) 



in all three cases. 



III. NUMERICAL RESULTS 
A. General remarks 

A detailed analysis of Eq. (21) reveals that for a sys- 
tem with a vanishing current the third line of the matrix 
equation vanishes identically. However, for these states 
the correct value of the current is already obtained at the 
level of SDFT as a natural symmetry constraint. In fact, 
using zero vector potential as the initial value, one can 
show that it remains zero at each iteration. Hence, one 
recovers the original SDFT result for non-current car- 
rying states. 13 On the other hand, for current-carrying 
states the xc vector potential is always non-vanishing 
even if one chooses a vanishing vector potential as the 
initial value. 

A closer inspection of the KLI equations shows that 
they become linearly dependent in the asymptotic region 
and therefore do not have a unique solution. In our nu- 
merical procedure, we take a pragmatic approach to the 
problem of linearly dependent KLI equations and add a 
very small positive constant to N(r) in Eq. (26). As the 
consequence, the limit becomes A xc {r) 0. In addi- 
tion, we impose v xc ^ a (r) —\Jr. This procedure also 
limits the possible appearance of numerical artifacts in 
the KLI potentials resulting from a finite basis-set. Such 
difficulties have also occurred for open-shell atoms. 13 ' 23 
Although we face similar problems in QD calculations 
(see below) , we have confirmed that the evaluation of the 
total energies, densities, and currents is not considerably 
affected. A further analysis is presented elsewhere. 24 In 
the context of the full solution of the OEP equations, 
problems in the computation of the effection potential 
due to the use of a finite basis-set have been recently an- 
alyzed in several works, and different possible solutions 
have been proposed. 25-29 



B. Examples 

Figure 1 shows the total energy of a six-electron QD 
(u>o = 5 meV) as a function of Bq- The kinks cor- 
respond to changes in the ground-state configuration 
(L Z ,S Z ). Apart from the fully-polarized (S z = 3) states, 
the EXX energies (dotted line) are considerably too large 
when compared with the accurate QMC results (dashed 
line). 8 EXX also leads to an erroneous occurrence of 
the (—5,2) ground-state at Bq = 1.5. ..2.0 T. However, 
adding the LSDA correlation 30 post-hoc to the EXX 
energies (EXX+cLSDA) yields the correct sequence of 
states as a function of Bq. This is a major improve- 
ment over the cLSDA-corrected Hartree-Fock calculation 
which docs not give the correct ground states for a similar 
system. 31 As expected, the corrections given by cLSDA 
are largest for the unpolarizcd state (0, 0) and smallest 
for the completely polarized states (—15, 3) and (—21; 3). 
This is due to the fact that the electron exchange has a 
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FIG. 1: (color online). Total ground-state energy (minus 
in a six-electron quantum dot as a func- 
tion of external magnetic field (SI units). The results have 
been calculated using the exact-exchange (EXX), EXX with 
LSDA correlation (EXX+cLSDA), and EXX with the cor- 
rected LSDA correlation (EXX+c||LSDA). The LSDA and 
quantum Monte Carlo (QMC) results 8 are shown for com- 
parison. The arrows mark the points where the ground- state 
configuration (L Z ,S Z ) changes. 

larger effect on the total energy in systems with a high 
number of same-spin electrons. 

Despite the improvement of EXX+cLSDA over the 
bare EXX. the result is not satisfactory in compari- 
son with QMC: Figure 1 shows that the energies of 
EXX+cLSDA are consistently too low by 1.0 — 1.5 meV. 
On the other hand, the agreement between QMC and 
the conventional LSDA (dash-dotted line) is very good. 
Hence, taking into account that the EXX is expected to 
capture the true exchange energy by a good accuracy (the 
only deviation arising from the missing correlation in the 
self-consistent solution), our result demonstrates the in- 
herent tendency of the LSDA to cancel out its respective 
errors in exchange and correlation. This well-known er- 
ror cancellation is lost when adding LSDA correlation 
to the EXX result. As expected, the performance of 
EXX+cLSDA with respect to QMC is at its best in the 
fully polarized regime (Bq > 5 T), where the exchange 
contribution in the total energy is relatively at largest. 

As a simple cure to the error in EXX+cLSDA, we ap- 
ply a type of self-interaction correction as first suggested 
by Stoll and co-workers. 32 The LSDA correlation energy 
can be improved by 

-EcTlLSDA = -EcLSDA 

- J d 2 r{n T (r)e c [n T) 0]+n;(r)e c [0,n;]}, (31) 

where e c [rt|,nj is the correlation energy per electron 
in the two-dimensional electron gas. 30 Therefore, in this 
approximation, denoted as EXX+c| jLSDA, the correla- 
tion energy between like-spin electrons is removed. We 
emphasize that this contribution is non-zero in the ex- 
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FIG. 2: (color online). Angular components of the paramag- 
netic and diamagnetic currents, j p g and jde, and their sum 
for the (—21,3) state at Bo — 11 T. The dashed line shows 
the exchange vector potential for the same configuration. 



act treatment and thus cannot be neglected. However, 
within the LSDA it contains mostly self-interaction en- 
ergy. Now, we find that EXX+c||LSDA (solid line) is 
very close to QMC, and actually performs better than 
the conventional LSDA. 

Figure 2 shows the paramagnetic current j p g and the 
diamagnetic current jd6{r) = n(r)Aog(r)/m* at Bq = 
11 T for the (—21, 3) state. The total current jg = j p g + 
jd$ changes sign at r ~ 350 a.u. due to the existence of a 
single vortex at the center of the QD. We find the vortex 
solution in agreement with both LSDA and numerically 
exact calculations. 33 

In Fig. 2 we also show the exchange vector potential 
A x . The small kink at r ~ 1100 a.u. is due to a basis-set 
problem described in Sec. Ill A. The maximum of \A X \ is 
located near the edge of the QD at r ~ 700 a.u. However, 
its relative magnitude with respect to the external vector 
potential Aq is largest at r ~ 150 a.u., where we find 
lAz/A) I ~ 0.1. Despite the considerable magnitude of 
A x , we find that its effect on physical quantities like the 
total energy, density, and current density is practically 
negligible. In the case presented in Fig. 2, for example, 
the difference between SDFT and CSDFT total energies 
is ~ 0.02%. In the context of the OEP method, the 
minor role of the xc vector potential has been observed for 
open-shell atoms, 13 molecules, 34 and extended systems. 16 
Earlier QD studies in the level of LSDA have also led to 
similar conclusions. 8 

Finally, we point out that, in principle, a given func- 
tional should be evaluated with KS orbitals obtained 
from self-consistent calculations and not in a post-hoc 
manner as we have done in this work. However, the vari- 
ational nature of DFT implies that if one evaluates the 
total energy with a density which slightly differs from the 
self-consistent density, the resulting change in the energy 
is of second order in the small deviation of the densities. 
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IV. SUMMARY 

We have applied the optimized effective potential 
method in current-spin-dcnsity functional theory to two- 
dimensional systems exposed to external magnetic fields. 
We have observed that the bare exact-exchange result 
(within the KLI approximation) is not sufficient in find- 
ing the correct ground-state sequence as a function of 
the magnetic field, although a considerable improvement 
over the Hartree-Fock results is found. Adding the cor- 
relation energy in the form of the standard local spin- 
density approximation yields excellent agreement of the 
ground-state energies with quantum Monte Carlo results, 
if the spurious self-interaction error is corrected. More- 
over, within the specified approximations, we found no 
considerable differences in total energies and densities 
when comparing the results obtained using a full-fledged 
current-spin-density functional theory and a spin-density 



functional scheme modified to include the coupling to the 
external vector potential. 

Acknowledgments 

The authors like to thank H. Saarikoski and A. Harju 
for providing us with their data published in Fig. 3 
of Ref. 8. We are grateful to Diemo Kodderitzsch for 
the valuable discussions on the properties of the effec- 
tive potentials and the issue of numerical artifacts in 
the optimized potentials. We gratefully acknowledge the 
financial support through the Deutsche Forschungsgc- 
mcinschaft, the EU's Sixth Framework Program through 
the Nanoquanta Network of Excellence (NMP4-CT-2004- 
500198), the Academy of Finland, and the Finnish 
Academy of Science and Letters through the Viljo, Yrjo 
and Kalle Vaisala, Foundation. 



1 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

2 W. Kohn and L. Sham, Phys. Rev. 140, A1133 (1965). 

3 U. von Barth and L. Hedin, J. Phys. C 5, 1629 (1972). 

4 G. Vignale and M. Rasolt, Phys. Rev. Lett. 59, 2360 
(1987). 

5 G. Vignale and M. Rasolt, Phys. Rev. B 37, 10685 (1988). 

6 P. Skudlarski and G. Vignale, Phys. Rev. B 48, 8547 
(1993). 

7 A. Wensauer and U. Rossler, Phys. Rev. B 69, 155302 
(2004). 

8 H. Saarikoski, E. Rasanen, S. Siljamaki, A. Harju, M. J. 
Puska, and R. M. Nieminen, Phys. Rev. B 67, 205327 
(2003). 

9 T. Grabo, T. Kreibich, S. Kurth, and E. K. U. Gross, in 
Strong Coulomb Correlations in Electronic Structure: Be- 
yond the Local Density Approximation, edited by V. Anisi- 
mov (Gordon & Breach, Tokyo, 2000), pp. 203-311. 

10 S. Kiimmel and L. Kronik, Rev. Mod. Phys. 80, 3 (2008). 

11 J. D. Talman and W. F. Shadwick, Phys. Rev. A 14, 36 
(1976). 

12 S. Sharma, J. K. Dewhurst, C. Ambrosch-Draxl, S. Kurth, 
N. Helbig, S. Pittalis, S. Shallcross, L. Nordstrom, and E. 
K. U. Gross, Phys. Rev. Lett. 98, 196405 (2007). 

13 S. Pittalis, S. Kurth, N. Helbig, and E. K. U. Gross, Phys. 
Rev. A 74, 062511 (2006). 

14 K. Bencheikh, J. Phys. A: Math. Gen. 36, 11929 (2003). 

15 S. Rohra and A. Gorling, Phys. Rev. Lett. 97, 013005 

(2006) . 

16 S. Sharma, S. Pittalis, S. Kurth, S. Shallcross, J. K. De- 
whurst, and E. K. U. Gross, Phys. Rev. B 76, 100401 (R) 

(2007) . 

17 S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002); L. P. Kouwenhoven, D. G. Austing, and S. 



Tarucha, Rep. Prog. Phys. 64, 701 (2001). 

18 E. Rasanen, A. Castro, and E. K. U. Gross, Phys. Rev. B, 
in print, cond-mat/0801.2373 

19 M. Ferconi and G. Vignale, Phys. Rev. B 50, 14722 (1994) 

20 S. Kiimmel and J. P. Perdew, Phys. Rev. Lett. 90, 043004 
(2003); Phys. Rev. B 68, 035103 (2003). 

21 R. T. Sharp and G. K. Horton, Phys. Rev. 90, 317 (1953). 

22 J. B. Krieger, Y. Li, and G.J. Iafrate, Phys. Rev. A 46, 
5453 (1992). 

23 D. Kodderitzsch (private communication). 

24 S. Pittalis, PhD Thesis, in preparation, Freie Universitat 
Berlin (2008). 

25 T. H. Burgess, F. A. Bulat, W. Yang, Phys. Rev. Lett. 98, 
256401 (2007). 

26 V. N. Staroverov and G. E. Scuseria, J. Chem. Phys. 124, 
141103 (2006). 

27 V. N. Staroverov and G. E. Scuseria, J. Chem. Phys. 125, 
081104 (2006). 

28 D. R. Rohr, O. V. Gritsenko, and E. J. Baerends, J. Mol. 
Structure: THEOCHEM 72, 762 (2006). 

29 A. Hesselmann, A. W. Gotz, F. Delia Sala, and A. Gorling, 
J. Chem. Phys. 127, 054102 (2007). 

30 C. Attaccalite, S. Moroni, P. Gori-Giorgi, and G. Bachelet, 
Phys. Rev. Lett. 88, 256601 (2002). 

31 W. Xin-Qiang, J. Phys.: Condens. Matter 12, 4207 (2000). 

32 H. Stoll, E. Golka, and H. Preuss, Theor. Chim. Acta 55, 
29 (1980). 

33 H. Saarikoski, A. Harju, M. J. Puska, and R. M. Nieminen, 
Phys. Rev. Lett. 93, 116802 (2004). 

34 A. M. Lee, N. C. Handy, and S. M. Colwell, J. Chem. Phys. 
103, 10095 (1995) 



